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Surface  disdrometer  measurements  were  also  simulated  accounting  for 
the  Poisson  and  Gamma  distributed  nature  of  the  drop  samples.  Sampling 
fluctuation  as  well  as  physical  variabilities  of  the  drop  size  distribution 
were  incorporated  into  disdrometer  simulations.  An  intercomparison  of  these 
simulated  variables  revealed  the  nature  and  reason  for  certain 
inconsistencies  reported  in  the  literature  by  various  authors.  Our  study 
showed  that  the  use  of  differential  reflectivity  does  indeed  improve  the 
characterization  of  the  median  drop  size,  and  that  its  use  in  rainfall  rate 
estimation  will  increase  the  accuracy  by  approximately  a  factor  of  2  over 
conventional  reflectivity  based  methods.  This  improvement  could  not  be 
detected  in  radar/disdrometer  intercomparison  experiments  because  of  the 
stochastic  nature  of  the  surface  data,  i.e.  ,  the  Poisson  and  Gamma 
statistics  of  the  raindrop  samples  introduces  random  fluctuations  very 
different  from  the  radar  fluctuations  thus  masking  any  potential  improvement 
predicted  from  theoretical  arguments. 

In  1987,  Colorado  State  was  chosen  as  a  Center  for  Excellence  in 
Geosciences  by  the  ARO  giving  added  impetus  to  multiparameter  radar  studies 
already  initiated  by  the  authors  of  this  report.  A  substantial  amount  of 
research  into  the  new  measureable  termed  differential  propagation  phase 
(^Dp)  and  *-ts  rai*ge  derivative  termed  specific  differential  phase  (K^p)  was 

initiated  using  time  series  data  from  the  CP-2  radar.  Using  experimental 
data  it  was  shown  that  in  rainfall  the  X-band  specific  attenuation  bears  a 
direct,  linear  relationship  to  S-band  ^R.  Simulations  showed  that  S-band 

4> Dp  could  be  used  to  predict  rainfall  attenuation  up  to  35  GHz.  S-band  4>Dp 

in  the  stratiform  ice  phase  up  to  1°  km  1  were  also  measured  and  interpreted 
as  being  due  to  oriented  ice  crystals  of  the  plate  type.  The  error 
structure  of  KDP  was  studied  using  simulation  techniques  developed  earlier 

and  it  was  shown  that  the  three  rainfall  estimators  based  on  reflectivity, 
ZDR  and  KDP’  termed  respectively,  R(Z)  ,  R(Z,ZDR)  and  P(KL)p)  havp  diff°*-er><- 

ranges  of  applicability.  For  rainrates  <  20  mm  hr’1,  R(Z)  appears  t:  have 
the  least  error;  for  20  <  R  <  70  mm  hr  1 ,  R(Z,Z^R)  appears  optimum  while  for 

R  >  70  mm  hr  ,  R(K^p)  is  optimum.  These  ranges  reflect  combined  error 

sources  due  to,  (i)  drop  size  distribution  fluctuations  and  (ii)  random 
errors  in  Z,  ZQR  and  K^. 

We  have  analyzed  the  axis  ratio  of  raindrops  as  function  of  their  size 
using  2D-PMS  probe  images  of  drops  in  convective  rainshafts.  Axis  ratios 
were  estimated  for  about  3500  drops  using  the  Fourier  descriptor  technique. 
Drops  with  diameters  less  than  4  mm  were  found  to  be  slightly  more  spherical 
than  would  be  expected  for  drops  in  equilibrium.  Oscillation  amplitudes 
were  found  to  be  typically  +10%  in  axis  ratio  for  light-to-moderate  rainfall 
rates,  and  such  oscillations  can  account  for  departures  from  equilibrium 
values.  We  have  determined  that  differential  reflectivity  can  be  used  to 
estimate  the  volume  median  diameters  of  almost-exponential  spectra  to  an 
accuracy  of  about  0.3  mm.  Without  this  information,  rainrates  based  on 
reflectivity  alone  could  be  in  serious  error. 
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